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The asymptotic wave function derived by Alt and Mukhamedzhanov [Phys. Rev. A 47, 
2004(1993)] and Mukhamedzhanov and Lieber [Phys. Rev. A 54, 3078(1996)] has been refined 
in the region where the pair (f3, 7) remains close to each other while the third particle a is far away 
from them (p a — > 00, r a /p a — > 0). The improved wave function satisfies the Schrodinger equation 
up to the terms of order 0(l/p a 3 ), provides the leading asymptotic terms of the three-body scatter- 
ing wave function with Coulomb interactions and gives further insight into the continuum behavior 
of the three-charged-particle wave function and helps to obtain 3 — > 3 scattered wave. This opens 
' up further ways of solving and analysing the three-body Schrodinger equation by numerical means. 
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• I. INTRODUCTION 

m : 

l/-) The quantum dynamics of three charged particles is described by Schrodinger's equation which should be supple- 

mented by proper boundary conditions. Merkuriev and Faddeev |l| claimed that the solution of this equation exists 
t-H ■ and is unique if the boundary conditions are known in all asymptotic regions. There are two types of the three-body 
scattering wave functions. The first type evolves from an initial three-body incident wave describing three incident 
particles in continuum. The second type of the three-body scattering wave function evolves from a two-body incident 
I wave corresponding to collision of a two-body bound state with a third particle. 

The three-body incident wave represents the leading asymptotic terms of the total three-body scattering wave 
function The knowledge of the three-body incident wave is important for many reasons: as a leading term 

I/"") , of the three-body wave function it can be used in calculations of the breakup matrix elements if the kinematics is 
such that the asymptotic region gives the leading contribution; the knowledge of the three-body incident wave is 
necessary in direct solution of the three-body Schrodinger equation for the scattering wave function of the first type. 
The asymptotic behavior of the three-body incident wave depends on the asymptotic region under consideration. In 
the asymptotic region where two particles, for example (3 and 7, are close to each other and far away from the third 
particle a, the three-body incident wave can be written as an asymptotic series in powers 1/ p a , where p a is the distance 
between the cm. of the system (/3,7) and the third particle a. The leading asymptotic terms 0(1) and 0(1/ p a ) of 
the three-body incident wave for charged particles have been obtained analytically in 0, Q . These asymptotic terms 
satisfy the Schrodinger equation up to 0(l/p 2 ). In this work we will derive all the leading asymptotic terms of the 
three charged particles incident wave of order 0(l/p 2 ). Combined with the previously derived terms of order 0(1) 
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and 0(1/ p a ) [2,11], they provide the asymptotic solution of the Schrodinger equation for three charged particles in 
continuum up to terms 0(1/ p a ). The terms 0(1/ p 2 ) satisfy first order differential equations. It is worth mentioning 
that the terms 0(1/ p^) are the next order terms compared to the three-body scattered wave which is 0(1 /i? 5/2 ), 
where R is the hyperradius. This term, as well as the leading-order asymptotic terms of the three-body wave function 
of the second type have been given in Refs. 0,0. Practical ways of extracting the scattering and breakup amplitudes 
using these asymptotic wave functions have been presented in Refs. @, 0] ■ 

The paper is organized in the following way. In Section I we introduce the three-body nomenclature and give the 
statement of the problem. In Section II we recall some of the important relations relevant to two-body scattering. 
In Sections III-V we present asymptotic solutions of the three-body Schrodinger equation in all orders which can be 
obtained analytically with the asymptotic method. Finally, Section F concludes the paper. 



II. STATEMENT OF THE PROBLEM 



We consider a non-relativistic three-body problem for charged particles of mass m a and charge z a , a = 1,2,3 in 
the continuum state. We follow the notations used in Ref. Q. The Greek letters stand for constituent particles of 
the three-body system or for the pair of two other particles. For example, a labels the particle or the pair (3 + 7. 
Such a supplemental notation is customary in few-body physics. The following conventional notations for the two 
body quantities are also used: A a = Ap~f, where a 7^ (3 7^ 7. The Jacobi coordinates are determined as follows: r Q 
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is the relative coordinate between particles and 7, and k Q is its canonically conjugated momentum. a a — m P m ~< is 

their reduced mass, mp^ = mp + m 7 . Similarly, p a is the relative coordinate between the cm. of the pair 7) and 

3 

particle a, and q Q is its canonically conjugated relative momentum. M a — m a mp^/M, M — ^ m v is total mass of 

v=l 

the three-body system. There are three sets of Jacobi coordinates r„, p„, where v = a, 0, 7. We frequently need the 
relations between the coordinates, and conjugate momenta for a channel v = 0, 7 and the corresponding a— channel 
variables. They are given by the following relations 



Pv \ _ I M-m„ tva M a \ I Pa 



(1) 



k " / V ~ 6 " a Alt ~ M-m„ / V k » 



(2) 



where v = /3,7 and the antisymmetric symbol e Q „ = — e ya! , with e Q;y = 1 for (a,v) being a cyclic permutation of 
(1, 2, 3), and e aa — 0. The motion of the three particles is described by the Schrodinger equation in the configuration 
space 

{E - T ra - T Pa - V}^ a (r a ,p a ) = 0, (3) 

3 2 

where V = V v , V v = (r u ) + V^ (r u ). The Coulomb potential is given by V^(r a ) = 13 - ~" e , Z v e is the charge of 

particle v. Similarly, V„ is the nuclear potential between the particles of the v— pair, where v = a, 0, 7. T Ta = — qtt 3 ") 

is the kinetic energy operator for the relative motion of particles and 7, and T Pa = —^§f- is the kinetic energy 
operator for the relative motion of particle a and the center of mass of the pair (0, 7), respectively. 

Our aim is to derive the asymptotic behavior of three-body incident wave up to terms 0(1/ p^) in the asymptotic 
region £l a , where r a j p a —> and p a — > 00. This incident wave provides the leading asymptotic terms of the three-body 
scattering wave function of the first type. General asymptotic behavior of the three-body scattering wave function is 
given by [l| 



Pv 
v—ol, p. 7 

, -^3^3 AkK-% X ln(2« R) (a\ 
+ R 5/2 ■ W 

Here the first term is the incident three-body wave, the sum over v provides the two-body outgoing scattered waves 
and corresponds to the 3^2 processes. The last term descibes the outgoing three-body scattered wave. Also, 
r\ a = (Zp + Z 7 ) Z a e 2 M a /q a is the Coulomb parameter for the Coulomb interaction between particles a and the 
center of mass of the system + 7; the Coulomb pasrameter Ao is determined in Ref. Q. We use the system of units 
such that h = c = 1. Formally we can determine the incident wave as an asymptotic difference 



<Pu(r v )— 



Pv 



M 



3^3 ^inR-iXo ln(2« R) 



i? 5 /2 



(5) 



From this equation it is clear that the three-body incident wave is a part of the full wave function, which does not 
contain the outgoing two- and three-body scattered waves. For better understanding of the three-body incident wave 
we consider first the two-body case. 



III. ASYMPTOTIC TWO-BODY SCATTERING WAVE FUNCTION 



We will be referring to the two-body Coulomb scattering throughout this work. Therefore we present here some 
important relations for the two-body scattering. Let us consider two charged particles with mass rrii and charge 
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Zi e, i = 1,2, interacting via the pure Coulomb potential V = Z\ Z2 e 2 /r. Scattering of two particles is described by 
the Schrodinger equation 

{E-H}^ + \r)=0, (6) 

where r\ = Zi Z2 e 2 /i/fc is the Coulomb parameter, E = fc 2 /(2/i) is the relative kinetic energy of the interacting 
particles 1 and 2, H = — A r /(2^t) + V is two body Hamiltonian, and /i = mi m2/{mi + JTI2) is reduced mass of 
particles 1 and 2. For the pure Coulomb interaction case Eq.® can be solved analytically. Substituting 

^ +) (r)=e ik - r AT 1 F 1 Hr 7 ,l,iC), (7) 
into Eq.© gives the differential equation for the confluent hypergeometric function (also called the Kummer function) 

A r ik ■ Vr 



2/1 /! 



V] 1 F 1 (-ir),l,iO=0, (8) 



N = e 7r? '/ 2 r(l + zjj) is the normalization factor, and parabolic coordinate C, = ki — k • r. iFi(—iri, 1, i() is called the 
Kummer function because Eq. JSJ rewritten in terms of z = i £ becomes the Kummer differential equation 0, [j| : 

d 2 iFi(a,c;z) diFi(a,c;z) 
2 \-[c-z) - a x Fi(a,c;z) = 0. (9) 

Here a = — 277, c = 1. Note that the Kummer function iF\{a, c; is a regular solution at f = (or r = 0) of the 

Kummer equation. Correspondingly, ^^(r) given by Eq. (JJJ is the normalized regular solution of the two-body 
Coulomb scattering problem. The Kummer function can be expressed in terms of the Whittaker functions W\ ^(z) 
using Eqs. (9.220.3) and (9.233.2) [jjj: 



+ fp— ) e^Wi,) _L_ e ,C/2 W ^ irh0 (e-^i C). (10) 

Each term of Eq. IjlUI) also satisfies the Kummer differential equation © providing a singular solution. Substituting 
Eq. Ill Oil into the Kummer equation leads to the Whittaker differential equation for each term: 

Here A = ±(1/2 + i 77) and z = ±i £. Evidently that both Whittaker functions in Eq. (I10[l satisfy the same Whittaker 
equation because it is invariant under simultaneous transformation z — > —z, A — * —A. Coming back to the normalized 
regular solution of the Schrodinger equation we can present it as a sum of two singular solutions: 

^+ ) (r) = ^N 1 F 1 (-ir,,l,iQ 

= ^ 0) (r)+^ sc) (r). (12) 

The first singular solution, as we will see below, is the incident wave 

4 0) (r) = e ik - r ^(C), (13) 
and the second singular solution is the scattered wave 

^ C ) (r)=elk . r ^(2) (a (14) 

^ 1 )(C) = e¥(i C )-l e ^^ +i)0 (iC), (15) 
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Evidently that for r\ = the incident wave becomes the plane wave e lk r and the scattered wave just disappears. It 
follows from Eq. (9.227) [13 that (|z| > 0) 



die' 1 



t + z 



(17) 



and from Eq. (6.922.2) [T(| 



die' 1 



t + z 



(18) 



Taking into account the asymptotic behaviour at \z\ — > oo of the Whittaker function, Eq. (9.227) [l{ 

(A-l/2) 2 



Wx fl {z) 

we derive the asymptotic behavior of J-^(Q: 



z x e- z ' 2 



IChoo 



1 + W_ 



Correspondingly the asymptotic behavior of the Coulomb distorted incident wave 

1 

"it' 



^ 0) (r) r =°° e ik - r e^ ln C [1 + 0(-^)]. 



The asymptotic behavior of T^ 1 ' (£) is given by 



-vq In 2kr 



l + 0(-) 



where f c is the on-the-energy-shell Coulomb scattering amplitude: 



f° = -v 



r(i + irfi 



H) 



-vq nr}/2_ 



r(i - it,) 

The asymptotic behavior of the scattered wave is given by 

to 



j — iri In sin 4 

2fcsin 2 1 



,ifcr 

r 



''' e- 4 " ln2fcr [l + C(4)]- 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



Taking into account Eqs. I|12fl . I|21|) . I|24|l we get the asymptotic behavior of the Coulomb scattering wave function 
for a system of two particles in the coordinate space: 



V^ +) (r) r e * k 'V ln C [1 + 0(1)1 + fC^_ e - lV ln2kr M + (l)i. 

iQ r iC, 



(25) 



Note that this asymptotic behavior is valid only for |£| — > oo. For r — > oo it is valid for all directions in the 
configuration space except for the so-called singular direction, for which k • r = 1. Here a = a/a. 

One can observe a very interesting feature in the case of the two-body Coulomb scattering. The regular solution, 
according to Eq. (|12fl . consists of two singular solutions, incident and scattered wave, each of them also satisfies the 
Schrodinger equation. Correspondingly, the asymptotic Coulomb scattering wave function consists of two terms. The 
first one, e lk ' r e l,,ln ^ (1 + O(^) is the asymptotic form of e lk r J r ^ 1 ^(iC) and represents the Coulomb distorted incident 
wave. The Coulomb distortion not only generates a logarithmic phase factor r\ In £ as an additional phase factor to 
the plane wave phase factor k • r, but it also generates an infinite series in powers of This is in contrast to the 
two-body scattering problem for particles interacting via short-range potentials, where the incident wave is given just 
by the plane wave. The second term in Eq. H25J) is the asymptotic form for e lkr 'T^ («£). It generates the outgoing 
two-body spherical wave times the Coulomb scattering amplitude and also contains a factor which can be written as 
an asymptotic expansion in powers of l/£. 
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IV. ASYMPTOTIC THREE-BODY INCIDENT WAVE 



After consideration of the incident wave for the two-body case, it is easier to proceed to the incident wave for 
the three-body case. By definition, the incident three-body wave is the part of the total three-body scattering wave 
function of the first type, which does not contain two- and three-body scattered waves. We have shown that the 
two-body incident wave is a singular solution of the two-body Schrodinger equation . It is naturally to ask whether 
the three-body incident wave is a solution of the three-body Schrodinger equation . An educated guess tells us that 
the answer may be "yes" . First, the scattering wave function of the first type which consists of the three-body incident 
and scattered waves (two- and three-body) is a regular solution of the three-body Schrodinger equation . However 
there are also singular solutions of the three-body Schrodinger equation . And it is very plausible that the three-body 
incident wave is one of the singular solutions while the scattered wave represents another singular solution. However, 
we cannot prove it until an analytical expression for the three-body incident wave will be available. 

Our goal in this work is to derive the asymptotic incident three-body wave function in the leading orders 
0(1), 0(1/ p%), 0(l/p 2 a ) in the asymptotic region f2 a , where particles (3 and 7 are close to each other and far away 
from particle a. We will demonstrate that the terms of order 0(1/ p 2 a ) can be derived without explicit solution of the 
three-body Schrodinger equation . In principle the method we use can be applied to get even the higher order terms 
of the three-body incident wave but it is worth mentioning that the next order term in the asymptotic expansion 
of the three-body incident wave 0(1/ p^ a ) is inferior to the outgoing three-body scattered wave 0(1/ R b ' 2 ). Hence 
for practical applications one need to know the three-body scattered wave before getting the term 0(1/ p^ a ) in the 
asymptotic expansion of the three-body incident wave. A knowledge of the three-body incident wave up to terms 
0(1/ p^) allows us to write down the leading asymptotic terms of the three-body scattering wave function of the 
first type in the asymptotic region £l a up to terms 0(1/ p^)- Note that the expressions for the asymptotic incident 
three-body wave in two other asymptotic regions Qp and f2 7 can be derived by simple cyclic permutation of indexes 
a, [3 and 7. As we have mentioned earlier, the asymptotic incident three-body wave is the part of the total three-body 
scattering wave function of the first type, which does not contain two- and three-body scattered waves. This wave 
function should smoothly transform into the asymptotic incident three-body wave function in the asymptotic region 
ilo- This smooth matching is an important part of the boundary conditions that provides a unique solution. 

The leading asymptotic term of the three-body incident wave in f2o derived by Redmond jlj, |l2| is given by 

y ( S ) Jr a ,p a ) = e lU ^e 1 ^ ]J e^ 1 "^ (26 ) 

v=ot, 0, 7 

where 



Va = Z ? Z l e Ma (28) 

K a 

is the Coulomb parameter of particles (3 and 7, p a is the reduced mass of particles /3 and 7. It is the three-body 
Coulomb distorted plane wave. For practical applications Merkuriev 0, Garibotti and Miraglia extended the 
asymptotic Redmond's term by substituting the confluent hypergeometric functions for the exponential Coulomb 
distortion factors. This extended wave function, often called the 3C wave function, is given by 

*Sl +) (ra,Pa)=e ik - r «e i ^"« \{ F v (Cu), (29) 

v— a, /5, 7 

where 

F v (C)=N vl F 1 (-ir )v , l;tC)> (3°) 
iFi(~iT] u , is the confluent hypergeometric function and 

N u = e-* 1 »' 2 r(l + ir h ,). (31) 

Note that 



^ k£V (r Q ) = e 4k - r °^(C) 



(32) 
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is the Coulomb scattering wave function of particles (3 and 7 moving with the relative momentum k Q and is well- 
behaved even in the singular directions (d, < C for r u — + 00) where the Redmond's asymptotic term is not determined. 
If any of the particles is neutral, then the resulting asymptotic solution becomes the plane wave for the neutral par- 
ticle and the exact two-body scattering wave function for the charged pair. However, neither Redmond's asymptotic 

term ^^^(r a , p a ) nor the 3C wave function (*a,Pa) are asymptotic solutions of the Schrodinger equation 

in the asymptotic domains Vl v , v — a, j3, 7. Redmond's asymptotic term, by construction, satisfies the asymptotic 
Schrodinger equation up to terms 0(l/r a , 1/rjg, However, in the asymptotic region, Q u , the distance between 

the particles of pair v is limited: r v < C. Hence the terms 0(l/r v ) are not small and the potential V„ in the 
Schrodinger equation has to be compensated exactly rather than asymptotically as happens when we use the Red- 
mond's asymptotic wave function in Hq. In the 3C wave function two very important effects are absent. Consider, 
for example, the asymptotic region f2 Q . In this region r a << p a . Hence the two-body relative motion of particles 
(3 and 7 is distorted by the Coulomb field of the third particle a . The second evident defect in the 3C function 
is the absence of the nuclear interaction between particles (3 and 7 which can be close enough to each other in fl a . 
Nevertheless, the 3C wave function can be used as a starting point to derive the the leading asymptotic terms of the 
three-body incident wave in Q a 2, 3], because this asymptotic three-body incident wave should match the Redmond's 
asymptotic term in f2o- We will demonstrate now how important the condition of the matching of the asymptotic 
wave functions is on the border of different asymptotic regions [2J . 

In the Redmond's asymptotic incident wave three logarithmic phase factors appear, one phase factor for each pair 
rather than two phase factors in the factorized solution. It is a very important conclusion. In all the conventional 
approaches for breakup processes, including coupled channels codes like FRESCO, the three-body scattering wave 
function is approximated by the factorized one. From the consideration above, it is clear that if Coulomb interactions 
are important, such an approximation is not accurate. If the interactions are short-range, the factorized solution 
matches the asymptotic solution in Qq and is justified in the asymptotic region £l a . It has been shown in that 

the actual asymptotic solution of the asymptotic Schrodinger equation (r Q , p a ), which matches the Redmond's 

asymptotic term in Qq, cannot be written in a factorized form and has a quite complicated behavior. In 0, Q all 

the leading asymptotic terms up to 0(1/ p a ) of the asymptotic wave function ( r a,Pa) have been derived in the 

asymptotic region Q a . In this work we will present a derivation of the expansion of the asymptotic wave function, 
S q Q (r Q , p a ), up to terms 0(1/ p a ). The derived asymptotic expansion contains all the terms 0(1), 0(1/ p a ) and 
0(l/p 2 a ). Since we are looking for the terms 0(l/p a ), we need to keep the terms up to 0(l/p a ). Instead of the 
asymptotic expansion of the Coulomb potentials Vp(vp) and V^(r 7 ) in terms of l/p a , we will start our derivation 
from the exact three-body Schrodinger equation The terms of 0(1/ p a ) will be dropped later. The asymptotic 
wave function in fl a should match the asymptotic wave function in The 3C wave function satisfies Eq.Q up to 
terms 0(l/r a , l/p a ) and we can use it as the initial wave function. However, this wave function should be modified 
to satisfy the Schrodinger equation in f2 Q . Note that usually in the literature it is assumed that the Redmond's 
asymptotic term satisfies the Schrodinger equation in f2o in the leading order only. First we will show that the 3C 
wave function satisfies the Schrodinger equation in f2g up terms of order 0(l/r^). To this end we just substitute the 
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3C wave function l|29[l into the Schrodinger equation 

(E - T ra - T Pa - V)[e^+^^ ^ kQ (r Q )^(r^)^ k7 (r 7 )] 

ifa a -i E Vp« In J • Vp„ 



Here 



Taking into account 



we derive 



2M a M a 

^k Q (r Q ) 



2^g ^ 



A P0 T=a,7 



2Mp Mp 



Vk^(r/3) 



+ 



^Pt . LO=a,l3 

2M 7 M 7 



]^(r 7 ). (33) 



(r w ) = JV„F(-* fZl/ , l;iC v )=^ 1) (iO+^ 2) (*C/). (34) 



+ ik w - V r. _ K] {rv) = Q (35) 

2^ 



(E - T ra - T Pa - V )[e^+^ <^ kQ (r a )^(r^)^ k7 (r 7 )] 

= 0(l/r2,l/r3,l/r3). (36) 

We did not use any approximation to get Eq. 1361 Thus the 3C wave function indeed satisfies the Schrodinger equation 
in S7o up to the terms 0(\jr 2 al l/r%, 1/r 2 ,), i. e. after substitution of the 3C wave function into the Schrodinger equation 
all the terms of order O(l) and 0(\/r a ) are exactly compensated. Hence the 3C wave function can be used as a 
starting wave function with a proper modifications to look for an asymptotic solution in fi Q . Taking into account 
Eq. H34f> we can rewrite the 3C wave function as in a form which is suitable for consideration in the Cl a asymptotic 
domain: 

*£ C qi +) ( r «>P«) =e lk - r "e 4 " 
x [jf\i{f,)JfHity)N a F a (i&) + ^(K^HK^NaF^iCa) (37) 
+Tf\i^)T^(iC 1 )N a F a {i<: a ) + Tf\iCp)T\ 2 \i^)N a F a (iC a )]. 

Here, asymptotically, for — » oo, the first term F^(^„) ~ O(l) and the second term T^^) ~ 0{1/C, V ). Hence 
in the asymptotic domain fl a and v = (3, 7, can be treated asymptotically while ifk a (r a ) should 

be considered explicitly, because f2 Q includes the region, where r a is limited. Moreover, in the asymptotic region 
Q a the relative motion of particles (3 and 7 is distorted by the third particle a due to the long-range Coulomb 
interaction. It means that the wave function of the relative motion of particles (3 and 7 in fl a will be different from 



8 



the wave function e lka ' r ° N a F a (i£ a ) describing the relative motion of particles (3 and 7 in the absence of the third 
particle. Since interacting particles j3 and 7 can be close to each other in f2 a , their nuclear interaction should also 
be taken into account. Following Q we replace each N a F a (iQ a ) in Eq. by the corresponding unknown function 

tpa (raj, n, m = 1, 2: 

*S Q +) (r Q! p Q ) = e ' k - r « 

x[^ 1) (i^)^ 1) «7)^ U) (r Q , Pa) +^ 2) (*C/3)^ 1) «7)^ 21) (r Q , Pa) 
+^ 1) (zC /3 )^ 2) «7)^ 2) (ra, Pa) +^ 2) ( i C/ 3 )^ 2) « 7 )^a 22) (ra, Pa)]- (38) 

Derivation of </5 a rl " l ' l (r Q ), n,m — 1,2 is our final goal. Now we substitute Eq. (|38|l into the Schrodinger equation 
©. When substituting Eq. I|38|) into the Schrodinger equation we assume that each term of the sum (|38f) satisfies 
the Schrodinger equation . Moreover, as we will see, each function </? Q (r a , p a ) depends on the preceding functions 
J r p l \iCf3)J~'y n \i(~ ( ) where n,m — 1,2, i.e. for each term in (|38|l the modification is different. We also take into 
account that 

(^A r „ +i — k„-V r „ -V^F^XKv) =0. (39) 

£j£v fly 

Substitution of the first term of Eq. (|38f) into the Schrodinger equation generates the equation for (r a ): 

H Z/I a ZM a [l a M a 

— £ V r Q ln^(iC,)-V, a + ytY1 V p« ^i 1] {Ku)-V Pa - y Q (r Q ) + 

u=p,'y v—p^ 

[l a ivi a 

x^ n) (r a , p a ) = 0. (40) 

Since particles /3 and 7 are allowed to be close in fi Q their interaction potential is given by the sum of the Coulomb and 
nuclear potentials. Now we will simplify this equation by dropping all the terms 0(l/p Q ) and explicitly compensate 
all the terms 0(1), 0(l/p a ), 0(1/ p^). We consider only the nonsingular directions, i. e. k„r„ ^ 1, u = /3, 7. To 
analyze the fifth term in the brackets we use equations 

Fl X) (Ku) = H l) «u)[l + 0(1/®], (41) 

f{i) {iCv) = e iv^C % (42 ) 



mp 1 k v ri (1 - k„ • v v f 



(43) 



(44) 



Pa 1 — e a2/ k„ • /3 Q P a 



(45) 



Vr -Wv-i -2 



771^ 1 ^otvpa 



+ 0(-t). 



" m /?7 *V P 2 (1 - ■ p a ) 2 p; 



(46) 
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To estimate the sixth and the ninth terms we use equations 



V Pq ln^CiGO = if*, - e„« k ". T " + 0(i) (47) 
^ 1 - k„ • f „ 

1 k„ — £ a vPa .-./ 1 n //irA 

= ^^ Q ~ hO(— ). (48) 

Pa 1 - € Q „k„ • p a P a 

To estimate the eigth term we use equation 

™/3 7 r„ i _ k„ • r„ 

Note that in Q a radius r a is limited a priori (more strictly, it is allowed to grow but slower than p a ). That is why 
we cannot use an asymptotic expansion in terms of 1/Ca m the asymptotic region fl a . Eqs 1(441) . 1)45(1 ■ l(46() and 1481) 
are valid only in while Eqs l|4H[l. ((47(1 and 1(49(1 are valid both in fio and £1 Q . 

Thus we reduced a three-body problem in the asymptotic domain Q a to a two-body problem: we need to find a 
solution of Eq. ((40|) . which describes the relative motion of particles f3 and 7 in the presence of the third particle 
a, which is far away, but it still distorts the relative motion of particles /3 and 7 due to the long-range Coulomb 
interaction. This distortion results in the dependence of (p a (r a , p a ) on p a . When p a increases this distortion 
should be weakened. Hence, Lp a lx \v a ^ p a ) actually depends on l/p a and 

V Pa ip^(v a . iPa )^^. (50) 

P a 

Because of that we may drop the second and sixth terms in Eq. 140(1 and rewrite it in the form 
[^-A r „ +i — k a -V ra +i-J-q Q -V PQ + — Vr.ln^^tW-Vr. 

-V a (r a ) + — V r „ hijf\i{p)-V ra ln^ 1) «7) 

Ma 

+ lk Vp ° ^^(^^P- ln ^ 1) (*C7)]^ U) (ra, Pa) = 0. (51) 

The last two terms are of 0(l/p a ). Note that to satisfy this equation up to terms of 0(l/p a ) all the terms of 0(l/p a ) 
must be compensated. Taking into account Eqs ((44(1 and 1(46(1 we can rewrite Eq. ((51(1 as 

[— A r „ + i — k^ (p a )- V rQ + i — q« • V Pq 

+ — E ^( r «-VpJ(V PQ ln^(zC, Q )-V r J-y a (r Q ) 
+ (£0 a £f a — + 4r) V p« ln^ 1) «/3a)-V PQ ln^W« 7Q )] 

771/3 7 -<"a 

x^ 11 )(r Q ,p Q )-0(l/p3). (52) 
We introduced here a new local momentum 

k W=k a -iX: — I^V^lnjW^^+u^ k y -c a ,,p„ 

„~^ 7 m /*r ^ Pa (1 - e av K ■ p a ) 2 

Note that variables V rQ and V Pcl are mixed up only in the fourth term of Eq. il">2l) . We are looking for a solution in 
the form 



Ml. , il.ll . , , , X(Pa), . ^i(l)( rQ ' P a ) , . 

<Pa '( r a, Pa) = ^„(o)( r a, Pa) (1 + — ) + ^ ' ( 54 ) 

Pa P a 
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where (r„, p a ) is a solution of 



[^-A rcv +! ik( u »(p a ).V ro -^(r,,)]^)^, Pa) =0. (55) 



2/ia " /X 

x(Pa) ~ O(l) and is a solution of the first order differential equation 

. 1 _ x(Pa) 
M a p a 

= -(ep a e ia — + V Pq ln^ 1) (i^ a )-V PQ ln£«(iC 7Q ). (56) 
Finally ^^m^ ' P") ~ ^(•"•) ^ s a soni tion of the inhomogemeous equation 

[-^A ra + * -i- k^-Vr a V a (r a )] pg^, p a ) = -i A q a • V Pa *£g)(r a ) 

E ^(r«-V p JV pQ ln^ 1 )( i C, Q )-V rQ ^ 1 ( 1 ( | ) (r Q ,p Q ). (57) 

Note that all the equations l(55*)l . and l(57jl are "two-body" differential equations. On the left hand side they 
contain gradients and Laplacians over only one of the variables, r Q or p a . Therefore these equations can easily be 
solved numerically. 

Now we consider the second term of Eq. . It satisfies the equation 

^ 2 )(^)^ 1 )(iC 7 )[^A ra + -l-A PQ + *^k a .V rct +*-Lq Q .V Pct + 

— [V ra lnrf\i(p) + V ra lnjrW(iC 7 )]-V ra 

+ i [Vp ° + ln^ 1) «7)]'Vp Q - V a (r a ) + 

— V r „ ln4 2) (iC/j).V Pa In^iiQ + -L V Pa lnJ^ 2) (i^)-V Po lnj^f)^)] 
x^ 21 )(r Q ,p Q ) = 0(l/p3). (58) 

Here, in the nonsingular directions (k„ • r„ ^ 1, v ^ a) 

<* s ~ „„ IM e «. p + 0(1,,. (59) 

1(1 —ir) u ) Q u Q v 

Also, in the nonsingular directions for v ^ a 

V ra \nF¥\iC„) = iV Va C„ + 0{\/r u )=i^k v {K-v u ) + 0(l/r y ) (60) 

171/3 7 

= i — — k v (k„ - e av p a ) + 0(1/ p a ) (61) 

"1/3 7 

and 

V Po InT^(iC) = iV PQ C, + 0(l/r„) - ie„a(-*A + k„) + 0(l/r„) (62) 

= ifc,, (p a - e aI/ k„) + 0(l/p a ). (63) 

When deriving l|58|) we took into account that 

(^-A rv +t — K-V rv -V v c )rt 2 \iQ=0. (64) 

To get an asymptotic equation from Eq. I|58|l which is valid up to 0(1/ p a ), all the coefficients of 0(1), 0(1/ p a ) 
and 0(1/ p a ) should be kept in the left-hand-side of the equation. Since in the nonsingular directions in fl a region, 
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J-"jp(iC/3) ~ 0(1/ p a ) only coefficients of 0(1) and 0(X/p a ) in the brackets of Eq. I|58|l should be left. Taking into 
account Eqs JHJ), (JJH} and <E3) we get 

[-L A r „ + i — k^(p a )-V ra - V a (r a ) 
zp a p a 

1 771 /3 1 1 foil 

+« r~ fc /5 — r a - p Q (p a -r a ))- V ra + i — qL 21) -V Pa (65) 

Ma Pa M a 

-ie al3 — k B (kp - e Q/3 p a ) ■ V Pa lnJFM(iC 7a )] ^ 21) (r a , p a ) = 0(l/p 2 ). 

7Tl a 

Here V Pn In^y 1 (zC 7 «) is given by Eq. l|45f) . We also introduced new local momenta 

k( 21 >(p Q ) = k a + k B fa - e aB p a ) +i (ir, + 1) - k/3 ~ £ ?^ a , (66) 



and 

qi 21) (Pa) = q Q + fc/9 (p a - e a3 k 8 ). (67) 
We also took into account that for v =/= a =/= t, i> t, e UT , e ua = — 1, and 

-e a7 (k/3 - e a fl p a ) + — — (p"a - e a flka) = -e a g (kg - e a g p Q ). (68) 

ttiq 7 M a m a 

We are looking for a solution of Eq. I|66|l in the form 

(21) / \ 

^^Pa) = V ^hr aiPa ) + Wr "' PaJ , (69) 

Pa 



where ^cf(o)( ra) S£L tisfies 



Mo 

Finally V Q (i)( r 

) ~ O(l) is a solution of equation 

r 1 . 1 



[^A ro +J -l k ( 21 )(p Q ).V ra - ^(raJ^SfgjCra, p Q ) = 0. (70) 



-1 



A r „ +i— k^p^-V^ -7 a (r a )]fl(r ai p a ) 

1 mjg 

Pa (21) _ (21) , x 

^-q a j -v PQ v Q(0) (r a , p Q ) 



2 

[« — ^2— (r Q - p Q (Pa-ra))'Vr a ]^^)(rc p a ) 



+«£a/3 ^ (k/3 - £ Q/ 3 p a ) ■ V Pa lnJF^^^a)^ (o)( r a, Pa)- (71) 



m, 



Since in Eq. (|71|l we keep only terms of order 0(1/ p a ) local momentum k£ 21 ^(p a ) can be replaced by 

k a(0) = k a + — k f3 ( k /3 -t a Pa)- (72) 

A formal solution of Eq. (JJTJ is 

O'a, A.) = Pa) + e - k ^(->- / d^G(r Q , 

[-[« 77 T~ k ( r a - Pa (PaT a ))-V ra ]^ 2 (l(r^, p Q ) 

Ma TO /3 7 

-i£a/3 — ^ (k/j - e Q/3 p Q ) • V Pa ln^W^' )^ 2 f 1 l(r^, p a )], (73) 
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Here V^d) ( r a j Pa) is a solution of the homogeneous Eq. (|70l) . 

The third equation for tfa 2 \r a , p a ) is obtained by substituting the third term in to 0. Following the same 

steps, which we used to derive the second equation, or just interchanging (3 <-> 7 in (|58J) we find ipa 2 \r a , p a ) in the 
following form: 

- ,,(12) , „ s , Ml) 



i L "i.-. l .p„i = g ) (, a , Pa ) / a( ; a,Pa) l «7-h 



Pa 



where (r a , p a ) is a solution of 

[-^A ra + J — k£ 12 )(p Q )-V rQ -V a (r a )]£(r a , p a ) = 0. (75) 
We can derive a similar to Eq. 1|71[> equation for ^l(r a , p Q ) which has a formal solution 
^ 1 g ) (r QlPa )=^ ) (r Q) p Q ) + e -^ 1) ^)- r ° | dr^G(r Q ,Oe k 

\ Yf\ n 1 1 {'\ r )\ 1 

[-[i T~ k P ( r a~ Pa (p Q -r Q ))-V r J^ ,L(r a , p a ) 

Pa TO /3 7 

-i£ a p— kpQzp -e a(3 p Q ) • V Pa lnJ^ 1) (iC 7Q )^ 1 ( 2 ( ] ) (r Q , p a )]. (76) 

The fourth equation can be derived after substituting the last term of Eq. I|38|l into Eq. © and it is automatically 
satisfied up to the terms of order 0(l/p a 3 ) in £1 Q because the product (iC,p)J-^ (i£ 7 ) = 0(l/p a 2 ). The fourth 



term in Eq. (|38|l leads to an equation for ipa 2 ^ (r Q , p a ) 



+ — [V rQ \n.Tf\i<#) + V r „ ln^ 2 )(iC 7 )].V rQ 
[V Pq ln4 2) ^) + V Pq In^f(iC 7 )]-V Pa 



-V a (r„) + — V ra ln4 2) (iC^)-V ra ln^ 2 )(iC 7 ) (77) 

f-t>a 

+ W a Vp « ln ^ 2) ( l C/5)-V Pcv ln^f(iC 7 )] x ^ 22 )(r Q ,p Q ) = 0(1/ pi). 

Using the same arguments we have used before, we may drop all the terms containing derivatives over p a when looking 

(22) 

for a solution in leading order. Then the equation for ip a reduces to 

[-^A ra +z — k^(p a yV ra - V a (r a )}^ a 22 \r a , p a ) = 0, (78) 

with a local momentum 

k^ (p a ) = k« + V ^- k v (K -e av p a ). (79) 

If V a is a pure Coulomb potential, = V£ , then Eqs l|55|) . 17UI) . (I78|l have the following solution 

^■>(r , pa) = N^\ Pa )F{-i4^{ Pa WC {ij) {Pa)), (80) 
Here, i = 1,2; j = 1,2 and N a ( Pa ) is defined as 

iVW)(p«) = e-^ W, (^/a r (i + ^W)( Pa )) j (81) 
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where r,™ (p a ) = Wj^f, and C fo) ( P«) = k<£ j) ( p a )r a - (p a ) ■ r a . 

K a {Pa) 

If V a is not a pure Coulomb potential, then the differential equations above, which parametrically depend on p a , 
should be solved numerically. Since all equations are of the two-body type, numerical methods are well developed 
and have been in use for a long time. They can be applied to solve the differential equations above as well. All the 
solutions found this way are valid in all directions of the asymptotic region fl a except for singular directions. 

Thus, returning to Eq. I|38l) we can claim that, having derived all four wave functions ^J, (r„ , p a ), i, j = 1, 2, we 
know the asymptotic behavior of the three-body incident wave of the scattering wave function of the first type in the 
asymptotic region Q, a up to terms 0(l/p^). 



V. GENERALIZED ASYMPTOTIC SCATTERING WAVE FUNCTION VALID IN ALL REGIONS 

Q v , v = a, /?, 7 

Now we are in position to present a generalized asymptotic scattering wave function which satisfies the Schrodinger 
equation up to second order and which is valid in all the asymptotic regions: 



(a/3 7 )(+) 
k Q q Q 



After substituting i|82|) into J2J and dropping the higher order terms we get, 

{E T ra T p - Q - V}[^^ +i ^^Jr a )^(r^ (r 7 )] 



ik Q -r Q +iq Q -p c 



ykJ r «)<^L ( r 7)[ 



^r a , ^k a • \7r 

2p a 



+e ik ^ r " + ^-"^ a (r a )^(r 7 )[- r 

+e ,k T .r T+ , qT .p 7 ^ (ra) (r/3)[ : 



2pp 



2p, 

' 0(^r, jr, r a ,rg,r 7 G ft 

P 7 

0(pr,pr), r a ,r 7 e 



Ma 

^ ■ y r/3 

P 1 



^7]<Pk,( r 7) 



0( 



1 1 



), r,,,^ 6 f2 7 



where the local momentum is given by 



t— a,/3,7 



(82) 



(83) 



(84) 



(85) 



In the asymptotic region f2oi each local momentum, k„, can be replaced by the corresponding asymptotic momentum, 
k„. In the asymptotic region Q a , Eq.JSJJ) reduces to the one quasi-two-particle differential equations: 



j A rQ ^ ik Q • \7r 



Mo 



V a ]<Pk (r ) = 0(— ,— ). 

773 r 7 



(86) 



Solution of this equation is evident and provides the Coulomb-nuclear scatttering wave function with the local mo- 
mentum k Q . Similarly we can get the asymptotic solution in leading order in the other two asymptotic regions fig 
and fL 



VI. CONCLUSION 



We derived the three-body asymptotic incident wave, which satisfies the Schrodinger equation in the asymptotic 
region f2„, v = a, (3, 7 up to terms of order 1/ p\- This asymptotic incident wave gives the leading asymptotic terms 
of the three-body scattering wave function of the first type and is an extention of the asymptotic wave function derived 
in 0,GJ. Equivalently, similar wave functions satisfy the Schrodinger equation up to 0(l/p^), v = j3, 7. It is worth 
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mentioning that the asymptotic solution satisfying the Schrodinger equation , in the asymptotic region fl v up to the 
0(1/ pi) can be found analytically - To find an asymptotic solution satisfying the Schrodinger equation in VL U 
up to terms of 0(1/ pf,) we need to solve two-body type differential equations numerically. The next order term in 
the asymptotic three-body scattering wave function represents the outgoing 3 particles — > 3 particles scattered wave 
and has been given in 0. 

The resulting asymptotic solution provides extended boundary conditions in all the asymptotic regions and can be 
used in the direct numerical solution of the Schrodinger equation or in approximate perturbation calculations as a 
leading asymptotic term of the three-body scattering wave function. 
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